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Abstract. In this work we provide an Aleksandrov-Bakelman-Pucci type estimate for a cer- 
tain class of fully nonlinear elliptic integro-differential equations, the proof of which relies on 
■ an appropriate generalization of the convex envelope to a nonlocal, fractional-order setting and 

, on the use of Riesz potentials to interpret second derivatives as fractional order operators. This 

» I ■ result applies to a family of equations involving some nondegenerate kernels and as a conse- 

^ j quence provides some new regularity results for previously untreated equations. Furthermore, 

, this result also gives a new comparison theorem for viscosity solutions of such equations which 

only depends on the L°° and L" norms of the right hand side, in contrast to previous com- 
parison results which utilize the continuity of the right hand side for their conclusions. These 
results appear to be new even for the linear case of the relevant equations. 



1. Introduction 

We begin this work with a very basic question to motivate our results. Suppose that L is a 
uniformly elliptic operator (L could be a second order or integro-differential operator) and that 
Uk are appropriate weak solutions (read viscosity solutions) of the equations 

\L{uk,x)=fkix) in Bi 

jn = onM"\Bi, 

with the additional assumption that < < 1. Then we pose the following question: 

Question 1.1. Under what conditions will it be true that \{x : fk{x) > 0}[ —)■ as /c —)> oo also 
implies that \\uk\\L°° — 0? 

(We will assume L is 1-homogeneous, and hence the constant function is a supersolution, and 
so always Ufc < 0.) 

In the case that L is a second order, uniformly elliptic operator, 

L{u,x) = aij{x)u^^x.{x), 
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(for Aid < (ttij) < Aid) the answer to the above question is indeed affirmative and is given by 
the celebrated Aleksandrov-Bakelman-Pucci (ABP) estimate which states that 

C(n) 

-mf{uk} < ||/fc||L" 
Bi A 

(where the equation is set in Bi C M"). 

The current understanding of this question is dramatically different, and not many results 
are known when L is a uniformly elliptic integro-differential operator, namely 



L{u,x)= {u{x + y)-u{x)-Du{x)-yl\y\<i{y))K{x,y)dy, (1.1) 

where A < K{x, y) |y |"+'' < A V x, y G M*". 

In fact, using current results the only occasions in which something could be said about the 
possibility of ||tifc||Loo — > would be when an explicit Green's function for the operator L in Bi 
is known to exist. For many applications and also the possibility of treating nonlinear equations, 
this is an unsatisfactory answer to Question ll.il 

In this work, we present a new Aleksandrov-Bakelman-Pucci (ABP for short) type esti- 
mate for subsolutions and supersolutions of integro-differential equations for particular kernels 
K(x,y). Namely, the kernels for which the result applies include (see Remark II. 5 1) 

K{x,y) = {2-a)f^^ for ctG(0,2), 

\y\ 

where A(x) satisfies a nondegenerate ellipticity condition 

Tt{A{x)) > a V X (1.2) 

(and is only measurable in x). Since these kernels are symmetric in y, we may rewrite L{u,x) 
as 

L{u,x) = {2-a) [ 6u{x,y)y-^^dy, a e (0,2) (1.3) 
JR" \y\ 

where we use the notation for second differences as 

6u{x, y) := u{x + y) + u{x — y) — 2u{x). 

We want to point out here the much richer structure inherent in elliptic integro-differential 
equations regarding uniformly and non-uniformly elliptic operators. A distinction must be 
made between the properties of the kernels which are used to construct operators such as (jl.ip 
and the properties of the operators, L, themselves. Indeed, it is very reasonable to consider L 
to be uniformly elliptic even if the corresponding K are not pointwise comparable uniformly 
to the standard kernel, (see [32] for recent regularity results in this direction). Such 

distinctions are by no means well resolved in the current literature, and this should be considered 
as evidence to the vast amount of work which is still to be done in this direction of integro- 
differential equations. 

Moreover, our result also covers fully nonlinear operators, such as 

F(u, x) = sup inf {La/3 (u, x)} (1-4) 
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where each Lq,^ is as in ()1.3p . with a fixed a for aU of the Lq/3. It gives new results for the 
Dirichlet problem: 

|F(n,x) = /(x) \i\D 

\u = g on M" ^ ' ' 

where F is given by (jl.4p . g is continuous, and D is an open, bounded domain. 

The ABP estimate is not only linked to convergence questions as explained above, but it 
also plays an important role in the regularity theory for solutions of equations related to (jl.3p - 
(jl.Sp . In particular, Theorem 11.31 below, implies both new comparison and regularity results 
for solutions of (II. 5p and gives a new proof of the recent Holder regularity results of [181 Section 
12] for the operators which are in the class exemplified by (II. 3p . Furthermore, regularity for 
equations with operators whose kernels obey the nondegenerate ellitpticity requirement in (|1.2p 
have not been previously attained, and so in Section [10] we briefly explain how to obtain Holder 
regularity results as a consequence of Theorem 11.31 

The paradigm of ellipticity taken in [18[ Section 3] (and [21[ Chapter 2 and 5] for second order 
theory) is that F in ([1.5p is considered elliptic if there are minimal and maximal operators, M~ 
and M^, such that 

M~{u -v)< F(u, x) - F{v, x) < M^{u - v). 

In this work, we use a more restricted version of these minimal/maximal operators than those 
in |18[ Section 3], which is reflected in the list of kernels covered as ([1.3p being much smaller 
than, but not exactly contained within those considered in [18] . Our extremal operators are 
defined as 

M-{u,x)= inf \{2-a) I 6u{x,y)^^-^^dy\ (1.6) 

A<Tr(A) and yl<AId \^ Vr- ' |y|"+'^+2 ^ ^ 

and 

M+(u,rE)= sup i{2-a) [ du{x,y)fL^dy] . (1.7) 

A<Tr(A) and A<Ald I JR" \y\ ) 

Supersolutions (respectively subsolutions) to (jl.5p for these F in the ellipticity family are char- 
acterized by the fact that they are automatically supersolutions of a minimal (respectively 
maximal) equation. Accordingly, we assume throughout this note u to be a viscosity superso- 
lution of a minimal equation with variable right hand side: 

{M-{u,x)<f{x) inBi 

\u>0 on M" ^ ' ' 

and without loss of generality we assume that / > and that u < in Bi. 

Remark 1.2. It is very important to remark that although the restriction A < Aid is necessary 
for the definitions of M" and M"*", the ABP-type result presented in Theorem 11.31 only depends 
on the value of A, and not that of A. This is to be expected for such a result, as seen Section [H] 
with the use of Lemma [ 



The main contribution of this work is to provide estimates on the infimum of u in terms of 
measure theoretic quantities of /, in particular for (jl.8p the and L" norms of /. The main 
result of this work is: 
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Theorem 1.3. Assume that u G P| L5C(]R") solves in the viscosity sense, that 

M~ is given in lil.6]) . and that f € C{Bi). Then there exists a constant, C{n), such that 

where C Bi is the coincidence set between u and a special envelope of u, made precise in 
SectionlM [34\ ). 

A few immediate remarks related to Theorem 11.31 are in order: 

Remark 1.4. For the definitions and basic properties of viscosity solutions of p.Sp and (|1.8p . 
the reader should consult the works: [3], [5], and \18\ Sections 1-5]. We emphasize the sign 
convention for subsolutions and supersolutions in this work corresponds to that of |18j and [21] . 
A more detailed history and presentation pertaining to viscosity solutions of first and second 
order equations can be found [25] . 

Remark 1.5. We would like to at least make an attempt to give a concise explanation as to 
the need for the peculiar restriction that A'I~ is of the form (jl.6p and hence why we can only 
treat kernels as in (jl.Sp . Roughly speaking, our approach starts with the function u, which 
solves a cr-order equation with a < 2, and takes its Riesz potential to invert the order by an 
amount 2 — a. Then to the potential, say P, of our original u, we can apply familiar second 
order techniques. In particular, we use the formula for the determinant of D^P(x) for those x 
with D^P{x) > 0: 

(det(L>2p(x)))^^" = -inf{Tt(AD^P(x)) : ^ > and det(A) = 1}. 
n 

When this computation is translated back to the original u, the derivatives on P can actually be 
transfered to derivatives on the kernel used to construct P, which results in a cr-order operator 
containing the term 

La{u, x) = {2-a) / Su{x, y) ,^,^ dy. 
JR" \y\ 

All of these steps are implemented carefully in Sections [3] - [5l 

Remark 1.6. Due to the slight degeneracy allowed by Tr{A) > A in the kernels, y"^ Ay |y| """^^^j 
comprising the definition of M~, Theorem ll.3l in fact leads to new regularity results. Specifically 
for functions, say w, satisfying boundedness of the maximal/minimal operators: 

M''{w,x) > -C and M-{w,x) < C, 

Theorem 11.31 provides a key step for the regularity theory, and this will be discussed further in 
Section [lOl 

The analysis of fractional order integro-differential equations has gained much attention lately 
(see Section [2]). It seems surprising, however, that despite so much interest, such estimates as 
those in Theorem 11.31 have not been obtained thus far. At least one need for an ABP type 
estimate for nonlocal equations was noted in Remark 3.3] regarding the ability to prove 
homogenization for stationary ergodic families of integro-differential equations via the methods 
of [22]. This stochastic homogenization result will be presented in [U]. Other potential appli- 
cations will be briefly discussed in Section [TTl In [181 Section 8] an ABP type result was proved. 
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but it involves the maximum of / over a finite collection of cubes, representing a Riemann Sum 
approximation to the usual ABP of second order theory. The result in Section 8] was 
sufficient for the purposes of regularity theory (for very general kernels), nevertheless it has not 
been sufficient to answer Question ll.il which is needed for the stochastic homogenization. 

To prove Theorem 11.31 we introduce new machinery, specifically a nonlocal and fractional 
order replacement for the convex envelope and the Monge- Ampere operator of a Riesz potential 
of u, which can be expressed as a nonlinear integro-differential operator in u itself. This auxiliary 
operator makes up for current lack of a definitive analogue of the Monge-Ampere operator 
(which has both a divergence and a non-divergence structure), which has been and continues to 
be a significant obstacle to understanding the more geometric aspects of the regularity theory 
for nonlinear integro-differential equations. 

We conclude this introduction with a brief outline of the sections of the paper. The two 
new tools mentioned in the previous paragraph, a replacement for the convex envelope and the 
auxiliary Monge Ampere operator, are explained in Sections [3] and [5] respectively. This new 
cr-order envelope solves a ci-order equation, just as in the second order setting, [211 Chapter 
3], with the convex envelope and a second order equation. However, in order to gain access 
to familiar geometric arguments involving the usual convex envelope, we must push the order 
of the envelope and its equation up to 2. This is done by taking the Riesz potential of the 
envelope, and its properties are presented in Section HI In Section we develop the Monge- 
Ampere operator of the potential as a cj-order operator acting on the envelope. In Section [H] 
we provide the final details of Theorem 11.31 In Sections [71 [H [9l and \W\ we present discussions 
involving respectively: the auxiliary operator of Monge-Ampere type, limits of these results as 
(7 — > 2, other useful theorems related to Theorem 11.31 including comparison, and a proof of the 
estimates of elliptic regularity theory which nearly identically follows the classical one of [2H 
Lemmas 4.5 and 4.6], from which the Harnack inequality follows easily. 

2. Background and Main Ideas 

2.1. Historical Background. Analysis of integro-differential equations is by no means a new 
field. It is intimately linked, via their infinitesimal generators, with modeling involving Levy 
or Levy-Ito processes- which are much richer and more general than diffusion processes giving 
rise to second order non-divergence equations. A key feature intrinsic in modeling with Levy- 
Ito processes and integro-differential equations is that they allow for long range interactions 
of various forms, and that the underlying stochastic processes can have jumps- in contrast 
to diffusion processes, which are continuous. There are many recent applications of these 
processes and their generators, and we list only a characteristic few: for particle systems and 
their hydrodynamic limit [27] . for financial modeling |23j and [23], for optimal control (also 
related to financial modeling) [39], [l5], and [H], for image processing [29] . 

There has been recently a growing interest in the analysis of integro-differential equations. 
The renewed interest seems not only due to the importance in modeling (Levy processes have 
been around for a very long time) but to advances in both the probabilistic and partial differen- 
tial equation analysis for treatment of equations related to (|1.5|) and (|1.8|) . Here we list only a 
few references for recent advances, and suggest the interested reader to consult those references 
contained within, as a complete list would be impossible. On the analysis/PDE side, integro- 
differential equations were brought to the viscosity solutions framework in [3^, [l5], and later 



6 



N. Guillen and R. Schwab 



in [2] , and the comparison theory for these equations was recently improved and re-presented in 
[S] , with the Dirichlet problem being considered in [3] . Holder regularity issues were considered 
in |42j for "non-divergence" form equations (cf. \42\ Section 3.6] for a discussion of "divergence" 
versus "non-divergence" in this context), and a series of works making a uniform theory for both 
Holder and higher regularity of more general versions of (|1.5p was done in [16], [T7], and |18] . 
Also by PDE methods, parabolic regularity was obtained for divergence form equations in |15j . 
Regularity for equations related to (jl.Sp was investigated by different, but still PDE, methods 
in [3]. Probabilistic analysis was used to investigate many integro-differential equations in both 
"divergence" and "non-divergence" forms to obtain various important results including Holder 
regularity, Harnack inequalities, and other finer properties in [T], [B], [7], [S], [S], and [TU] . 

2.2. Another Simple Motivation For (jl.Sp and Theorem 11.31 Let us now setup another 
very simple question related to that posed in the Introduction. It arises in the study of obstacle 
problems, and was utilized in [22] for homogenization (see |41j for the same use in a nonlocal 
setting). Consider the solution of an elliptic equation, 

F{v,x) = inBi 

v = ^ onW\Bi ^ ' ' 

and the solution to the corresponding obstacle problem with the same operator, u, solving 



u = sup{0 : F{(j), x) > in Bi and < in all of M"}. (2.2) 

A reasonable question is: what is the difference between u and vl Naturally, since u is a 
subsolution of the equation for v and they share the same data in M" \ Bi , then at least we can 
conclude v > u. But what about the reverse inequality? Ellipticity (see \18\ Definition 3.1] and 
[2H Chapters 2 and 5]) tells us that 

M"(n -v,x) < F{u,x) - F{v,x). 

Two of the key properties of the solutions to obstacle problems is that u is actually a solution 
to F{u, x) = whenever u{x) ^ 0, and that u inherits the equation from when it does happen 
that u{x) = 0, which is F{u, x) < F(0, x) (here we mean the operator applied to the constant, 
0, function). Therefore, the relevant inequalities are 

M^(n -v,x) < F{u,x) - F{v,x) < l[^^Qy{x)F{0, x), 

which by the boundedness of F{0, x) <C becomes 

M~ {u -v,x) < Cl{„=o} {x) in Bi 

u-v = G onM"\Bi. ^ ' ' 

After considering an / which is a continuous approximation from above of l|„=o}(2^); one can 
clearly see the question at hand and the motivation for this work: How can (|2.3|) . and hence 
also (|1.8p . be treated in a way which only depends upon measure theoretic properties of / and 
not continuity properties of /? 

This leads us directly back to Question 11.11 The natural setting for the homogenization 
problem will involve a sequence of Uk and v as above with the additional information that 
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\{uk = 0}| ^ as A; — )■ oo (see [HI Section 3, specifically Lemma 3.5]). The goal in such a 
situation is to be able to conclude that 

ll^fc ~ ''^IIl"" — as A; — > oo. 

Before this work was completed, there seemed to be very little results, in fact the authors found 
none, which could handle a situation in which the right hand side of (11. Sp would converge in 
any sense other than in (the relevant results in both [3] and [18] strongly require / to 
be continuous and converge to in to conclude similar statements). At least within the 
restricted class of operators presented here, Theorem 11.31 answers this question, which was not 
possible with previous known results. 

2.3. Theorem 11.31 As a Proof of Concept For Broader Results. As hinted in the in- 
troduction, the full class of nonlinear nonlocal elliptic operators is much richer than simply 
those given in our definition of M~ . In [18] , there is a version of which covers many more 
operators than the one appearing in p.6p : it is given as 

M,-,(., X) = / (2 - .) (x^^^^ - A<^^(|#1) iy (2.4) 

Jm." V \y\ \y\ J 

= mi {{2-a) [ 6v{x,y)-^^j^dy] . (2.5) 

\<a(x,y)<A [ J^n \y\^^ J 

Hence the results of |18j correspond to a much larger family of equations than the one here. To 
go even further, as done in the linear case considered in [8], one could use this same definition 
with measures, n{x,dy), "comparable" in a less restrictive sense to \y\~^~'^ dy as opposed to 
only those measures with a density as above, n{x,dy) = a{x,y) \y\~"'~'^ dy, for a uniformly 
bounded from below and above. 

A reasonable guess for an ABP type theorem applying to (|1.8p in these classes above would 

be 

-inf(«)< ^11/11^., (2.6) 

for some p > n, depending on a. At the time [18] and [IQ] were completed, it was not known 
whether or not such a result should be true for the general nonlocal ellipticity class or even a 
restricted class such as the one considered here. Now Theorem 1 1 . 3 1 indicates that at least some 
form of ABP type result holds for a restricted class of operators. This gives hope that ABP 
type estimates, such as those in Theorem 1 1.31 for the general class of nonlocal elliptic equations 
may still be true. Furthermore, it opens the door to answering the question of whether or not 
(j2.6p is appropriate to expect for (jl.8p . The moral of the story is Theorem 11.31 indicates it is 
not that ABP type results are incompatible with the intrinsic properties of (|1.8p . but more 
importantly that the existing machinery is incompatible with (jl.Sp . This illuminates one of the 
main difficulties in analysis of nonlocal equations: there is no known general framework to take 
the place of the very important Monge- Ampere operator and the gradient mapping of convex 
functions in the second order theory. 

2.4. Failure of The Convex Envelope. For the convex envelope, say F, the information 
in the proof of the original second order ABP estimate is naturally encoded in the set {x : 
det{D'^T{x)) > 0}. This matches very well with second order equations because det(L)^r) = 
whenever F and u do not coincide, and on the set where they do coincide, F inherits the 
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supersolution property of u simply by comparison. In the nonlocal setting, the function w = 
|x|" — 1 will solve (|1.8p for a > a with a right hand side, /, which will still be bounded and 
continuous. However for a < a < 1, the convex envelope of w only coincides with w at one 
point, X = 0. Hence there would be a contradiction with the usual ABP estimate which would 
read: 

-infM<C||/||in(|^=r})=0. 

The problem here is exactly that such fractional order equations allow for much more drastic 
bending of supersolutions than is possible in a second order setting. Therefore, the convex 
envelope is not well suited for a measure theoretic estimate, such as Theorem 1 1.31 For purposes 
of studying the regularity of (II. 5p with even more general operators than the F appearing in 
p.4p . the convex envelope was sufficient and strongly used in [18l Section 8]. In contrast to 
|18| . in this work we must construct a different envelope which will be better suited to handle 
such a function as w, above. This will be the content of Section [3] 

2.5. Main Ideas and Sketch of The Proof. Here we give a brief sketch of how Theorem 
11.31 is proved. The main ideas are the same as in the second order theory (cf. |2H Chapter 3] 
or |28^ Chapter 9]), but the machinery and implementation are a bit more involved. 

Everything starts with an appropriate envelope of u from below, which we will denote as T. 
It must be such that 

inf{n} = inf{r}. (2.7) 
The first key feature is the existence of an operator, which we denote as V^j, such that 

(-infrf<c/ {V^(T,x))Pdx (2.8) 

for some p possibly depending on n and a. At this point for the second order theory, we 
would have F as the convex envelope of u and Da would be det(Z)^F)^/"', in which case the 
previous inequality is a consequence of the geometry of convex functions (what is known as 
"Aleksandrov's estimate"). The second key feature is that the operator must satisfy for a 
special class of F, 

Va{T,x) <CM-{T,x). (2.9) 
Finally, the third key feature is that 

Da{T,x) < whenever r{x) ^ u{x). (2.10) 

This is essential so that all the contribution of 2?cr(F) can be ignored except for the contact set 
between u and F. This way, information about M~(u,x) < f{x) in the viscosity sense can be 
carried over to F via comparison, and the other values of ^'^(F) will not pollute the integral 
in the estimate (|2.8p . Essentially F acts as a test function on u, and the defining feature of 
viscosity solutions is that at those points where F touches u from below, 

M-(F,x) </(x). (2.11) 

If all of (gZD, (gSD, (1212D, ([2TTU|) . and (pTT|) can be satisfied (which is very nontrivial), then 
Theorem 11.31 can be proved. 

It turns out to be quite difficult to simultaneously achieve all three of the key features, (j2.8p . 
(|2.9p , and (|2.10p . This delicate balance is what leads to Theorem 11.31 only being proved for a 
restricted class of operators, instead of the much more general class of |18] . 
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2.6. Notation. We list here some notation which will be used throughout this work. 

(1) The second difference operator: 5v{x, y) := v{x + y) + v{x — y) — 2v{x) 

(2) The n — 1 dimensional sphere S^~^ = dBi C M" and its surface area ujn 

(3) The complement of a set, = M" \ ^ 

(4) The following universal constant will appear often (when n > 2), 



A{n, a) = 7r° 



" r(^ 



r(f; 



It can be showed that A(n, a) ~ a as a — )■ 0"'" (see [Ml Chapter I, page 44]). 

(5) The Riesz Potential of order a, Ka{y) ■= A{n,a) ly]""^" 

(6) The Fractional Laplacian (when n > 2) 

-(-A)"/^^;(x) = "^"+""^^ A(n,2-a) / 6v{x,y)\y\-''~'' dy. 



(7) The one dimensional Fractional Laplacian in the direction r € S" ^ , (for some universal 
constant (A(l, 2 — a) which we will not need to specify here), 

-{-A)'^/%{x) := ^^^±p^A{l,2 - a) [ {6v{x, st)) Isr^"" ds 

Jr 

(8) The inf-convolution of a function v, V£{x) := miy{u{y) + (2e)~^ |a; — yp} (See [301 
Equations (14), (15)], [35], [SSI Appendix]) 

(9) V is C^'^ from above at x (respectively from below) [181 Definition 2.1] if there exists a 
radius r, a vector p, and a constant M such that for all \y\ < r 

v{x + y) — v{x) — p ■ X < M |y|^ 

(respectively v{x + y) — v{x) — p ■ x > M |y|^) 

(10) Contact Set between u and its envelope, F (defined in Definition 13. 6p . 

Ku := {x : u{x) = T{x)} 

(11) The convex envelope in Br, v^^ , for a function u > in M" \ Bji is 

v^^{x) = sup{/(x) : / affine and / < -y in -B_r} 

(12) The extremal operators and ellipticity constants for the family governed by L of ()1.3p 
are M" and M+ defined in ([r6|) and (fTTTl) . 



3. The Fractional Envelope 

This section is dedicated to constructing a new envelope (as a replacement for the convex 
envelope) for the supersolution, of (jl.Sp which will be essential to proving Theorem 11.31 The 
main idea is to imitate the most important features of the convex envelope as they pertain to 
the second order theory, cf. [211 Chapter 3] and [131 Proposition 2.12 and Appendix A] for the 
fully nonlinear version and [281 Chapter 9, Section 1] for the linear version. The goal is for the 
new fractional order envelope to be in a class of functions for which M~ is comparable to a 
nonlocal version of the Monge- Ampere operator, and also to cause this nonlocal Monge- Ampere 
operator to vanish whenever u and its envelope do not touch. For the sake of explanation, let 
F be the convex envelope of u and M~ to be the second order minimal Pucci operator (see 
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Chapter 2, Section 2]). For F, the two requirements just mentioned above correspond to the 
two facts implied by convexity and the envelope property: 



JJ efc < ( — ^ e/c I < f — A4~(F) j for the eigenvalues, e^, of 



and 

det(L'^F(x)) = if u{x) ^ T{x). 
The appropriate analogs to our new envelope appear subsequently as Lemmas 15.81 and I5.10[ 

In order to define the new envelope, we introduce a matrix-valued integro-differential opera- 
tor, which will play the (auxiliary) role of an "integro-differential Hessian" . 

Definition 3.1. Let n > 2, o" € (0, 2), and let v he a hounded function such that 

\5v{x,y)\ 

\y\ 



then we define 



h^^y^ ^) := ijl±ZzMn±^A{n, 2 -a) ^^<5^(x, y)dy (3.1) 

\y\ 



The use of this matrix is dictated by our approach based on Riesz potentials in Section HI and 
it will become more clear in Section [5l By this we mean the following, if v is smooth enough 
and K2-fT is the Riesz Potential (see section [2T6|) then a consequence of Lemma [5TT] is 

D'iv * K2-.) = h,{v, x) + ^^^^ Id 

Remark 3.2. It is immediate that if u touches v from above at x and they are both smooth 
enough then we have the matrix inequality ha{u,x) > ha{v,x). 

Remark 3.3. As a — > 2~ (and thus a ^ 0), the above identity gives us 

h^{v,x)^DMx) + ^ld 



Remark 3.4. The fully nonlinear integro-differential operators that our methods handle are 
exactly those that can be written in the form 

I{u,x) = mfsup{TT{A''\x)h„{u,x) : Tt{A''\x)) > A and A^^x) < Aid}, 
b 

where a and b can be taken over arbitrary index sets. 

Moving forward, the "a-order envelope" of u is defined as the solution of an obstacle problem 
which is a generalization of the obstacle problem satisfied by the convex envelope. First we 
define an operator which will be a fractional order replacement of the first eigenvalue of the 
Hessian operator used in the construction of the convex envelope for second order equations. 
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Definition 3.5. Let Xi{B) = min{e : e is an eigenvalue of B} be the first eigenvalue operator 
of a matrix. E„ is defined as 

E„{v,x) := Xi{h^{v,x)} = inf x) r) • r}. (3.2) 

res'"-! 

Now, using E(j, we can define our ^^a-order envelope" . 
Definition 3.6. 

T^{x) = sup{v{x) : Ea{v) > in B3, and v <u1b^ in M"} , (3.3) 
and the contact set between and u will be denoted as 

Ku:={x: rl{x)=u{x)}. (3.4) 



Remark 3.7. The usual convex envelope used in the second order theory has an analogous 
structure as the solution to an obstacle problem. The interested reader should consult [38] for 
an exposition. 

By remark [3.31 this new operator, E„, does not recover the smallest eigenvalue of the Hessian 
of u in the limit a ^ 2, instead (see Proposition 18. ip 

lim E^{u,x) = Xi{D^u{x)) + -^Au{x). 

In particular, as cr — ?> 2^ the envelope T^ converges to the solution of the upper obstacle problem 
for the operator above with u as the upper obstacle. The solution to this problem lies above 
the convex envelope of u but it does not agree with it. 

Remark 3.8. All of the results of this section hold (with small modifications) if instead we use 

E:{u,x):= inf {-(-A)>(x)} 

Here (— A)^n is the one dimensional fractional Laplacian of u defined in Section [2.61 This is a 
much more drastic notion of envelope (such an envelope would be below the one defined using 
E(j, and touch n on a much smaller domain, which can be problematic). In principle, it is 
perfectly tailored to handle much more general kernels: 

K(y) = "'j-^} where a(ry) = a(y) V r € M. 

y TL-\~(T 

However, it is not yet clear how one can go about proving bounds for this envelope in terms 
of a "convenient" integral quantity (i.e. one that can be controlled by and integral of M~{T)). 
In contrast, the envelope we use admits integral bounds granted by its compatibility with the 
Riesz potential. Whether this argument work for the more "drastic" envelope is not clear and 
perhaps an entirely different approach is need, this question will be addressed in future work. 

Remark 3.9. Some nonlinear one directional operators related to this alternative operator E* 
have also been considered in [11]. In particular, the one dimensional Fractional Laplacian in 
the direction of Vn gives rise to a natural integro-differential analogue of the Infinity Laplacian. 

For the remainder of this section, some properties of F^ will be collected for later use. Also, 
we will dispense with the notation F^ and instead simply use F except in special cases which 
will be appropriately noted. 
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Lemma 3.10 (Compact Support of T). T = in M" \ ^3. 

Proof of Lemma Vj.liA This lemma is an immediate consequence of the equation satisfied by T, 
which is treated as ()12.4p in the Appendix. 

□ 

As a solution to an obstacle problem, T attains regularity from two sources. A one sided 
regularity of u from above is transferred to F via the obstacle, and further regularity of F from 
below is attained via the structure of the operator, E^j. We record the result as it pertains to 
F in Proposition 13.111 but the result holds in more generality and is of independent interest- 
we have chosen to present these related results in the Appendix, Section [T2l 

Proposition 3.11 (F Regularizes From Below). Let CC B^. If u satisfies respectively 

(i) -(-A)""/^!! < C classically in Bi, 

(ii) u is C^'^ from above with a bound, C 

then there exists a Ci depending on fl, a, n, X, A ||ti||Loo and C such that for a.e. x € O, 
respectively 

(i) 

< /i(j(F, x) < Cild and \\^\\H^{n) 1^ Cii 

(ii) 

mx,y)\ 



\y\ 



-dy < Ci. 



Remark 3.12. In particular, Proposition 13.111 ensures that for almost every x, the a order 
operators can be evaluated classically on F. 

Remark 3.13. Proposition 13 . 1 II would in fact be easier and stronger if the obstacle problem for 
F were posed in all of instead of B^. In that case, F would be C^'^ from above whenever 
u was C^'^ form above at well, and would follow from a straightforward adaptation of |43t 
Proposition 3.10] to the case of a concave nonlinear operator. 

Proof of Proposition \3.11[ The proof is an almost immediate consequence of Proposition 112. 1[ 
We simply note that the fact that E^{T,x) > in S3 implies that hu{T,x) > in B-^. Hence 
/lo-(F) has all non-negative eigenvalues and so Tr(/io-(F)) bounds all of them individually. Hence 
because 

Tr(/i.(F)) = -{-Ar/^T, 

the result (i) follows immediately from Proposition 112.1"! Result (ii) follows immediately from 
Corollary nisi □ 

The following fact is a straightforward consequence of the definitions of M~ and , we state 
it as a lemma without proof. 

f \Sv(xj y) I 

Lemma 3.14 (Ordering of Ea^ and M~). For any v which satisfies / ^' < 00, 

Jm" \y\ 

Ea{v,x) > ^ E„{v,x) < ^^M~{v,x). (3.5) 

A 
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Exactly as with Lemma [3.14t the following is a straightforward consequence of the definitions, 
and so we again omit the proof. 

Lemma 3.15. If u satisfies classically — (— A)""/^!* < C in Bi, then for almost every x € .63, 
the operators on T are classically defined and satisfy for all A > 0: 

Jr" \y\ 

To wrap up this section, we prove one last property of the envelope, T. This is the stability 
of the contact sets with respect to an increasing family of obstacle functions, U£. 

Lemma 3.16. If u as e \0, then lim supine C K, where = {x : Ue{x) = Tf.{x)} and 
K = {x : u{x) = T{x)} and T^, T are as in Definition \3.(A 

Proof of Lemma \3.1(A The increasing property of Ue implies for £2 < Ei that T^^ is an admissible 
subsolution below , and hence because T^^ is the largest such subsolution. Teg > . Thus 
Fe is also increasing, F. Now suppose x is a point in the set lim supj_^Q Kg . We can 

extract a subsequence indexed by £k such that n^^, (x) = F^^, (x). 

It must be shown that F(x) = u{x). We already know that F(x) < u{x) by definition of F, 
so therefore we will show F(x) — u{x) > 0: 

> F(x) -u{x) 

> T{x) - Ts^{x) + Te^{x) - Ue^{x) + Ue^{x) - u{x) 
>0 + + Ue^{x)-u{x). 

Letting A; — > cx) (e^ — )■ 0) yields the result. □ 

4. The Fractional-Order Potential 

In the previous section we defined and presented some properties of an important envelope 
which will be used in the proof of Theorem 11.31 However, there are still many difficulties to be 
overcome in proving an ABP type estimate for (11.80 . specifically one must relate the infimum 
of F to an integral of a quantity which is comparable to M~(F). In the second order setting 
(using a convex F), this is resolved by the very special geometry of the convex envelope which 
gives 

-mf{T}<C{n) [ det{D^T{x))dx. 

A key feature in the proof is that we can transfer this argument to the Riesz potential of F, 
which will solve a second order equation. Moreover, the usual argument will result in looking at 
the gradient measure, \VP{B-^)\, which is a nonlocal quantity of F. Inspiration for this choice 
came from the use of the Riesz potential in the setting of a nonlocal Porous Medium Equation 
treated in [19]. This method is also reminiscent of the regularity theory for the obstacle problem 
for the fractional Laplacian 03] , which starts with the observation that the fractional Laplacian 
of the solution (of order cr) solved an equation of order 2 — cr. 

As mentioned, we take our fractional-order potential to be the Riesz potential of the envelope 
(cf. [ini Chapter 5, Section 1] or [3ll Ch. 1, Section 1] ), which formally says 



p = (-A)-(2— )/2r. 
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Precisely we mean for this to be the convolution with the Riesz kernel, Ka '■= A{n, a) |-| 

P = T*K^, (4.1) 
for a = (2 — 0"), and the constant A{n, a) which is listed explicitly in section [2^ Hence 

(-A)P = {-^yl^T and (-A)(2-'^)/2p = r. 

The whole motivation of using the potential, P, of order (2 — a) /2 is to make sure that the 
usual 2nd order Monge- Ampere operator, dei{D^P), becomes a u-order operator when viewed 
as an operation on V. This way, there is hope that d.ei{D^P) will be comparable to M~r, 
and familiar techniques from the proof of the second order ABP estimate can be used (see 
Lemma 3.4] and/or [281 Lemma 9.2]). Indeed, the form of dei{D'^P) is investigated in Section[5] 
and will be shown to be a cj-order operator comparable to M~T (this new operator will not be 
an elliptic <T-order operator in general, but we will be working with a convenient special case). 
However, some additional difficulties are introduced by using P instead of F itself. It will be 
necessary to compare the infimum of P to the infimum of T in a uniform fashion which does not 
depend on the continuity of F or u. This is not to be expected simply from the Riesz potential 
itself, but must be deduced from the equation satisfied by F. The relationship between P and 
F is the main result of this section, which is presented in Proposition 14.71 

To begin the results of this section, it will be helpful to remark upon the regularity of P 
pertaining to that of F. These results are presented in [Hj, and so we do not provide a proof. 
They all follow from either of the two equations involving P: 

(_A)(2— )/2p = r 

or 

(-A)P = (-A)'^/2p_ 

Lemma 4.1 (Proposition 2.9 of |43|). The regularity on P inherited from F is 

i) whenever (2 — a) < 1, P G C*''" for any a < {2 — a) and 

||P||C0.C< <C(||P||oo + ||F||oo), 

ii) whenever (2 — o") > 1, P G C^'" for any a < {2 — a) — 1 and 

||P||cl,a < C(||P||oo + ||F||oo), 

and 

iii) P G C2'7 whenever F G C^^'^{W) (only for a <2). 

There are two more very useful consequences which follow simply from the equation satisfied 
by P. They will play a significant role later in both Proposition 14.71 and the proof of Theorem 
Ol 

Lemma 4.2. If xq G and I is a linear function, then P — I cannot attain a local minimum 
at xq. 



Proof of Lemma \4-'A noted above, P globally solves the equation 

AP = -(-A)^/2r_ 
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Furthermore, due to the fact that F satisfies the Dirichlet problem, (112. 4p . we know F(xo) = 
and without loss of generality F < in i?3, we also have that 

AP(x) = -{-Ay/^T{x) < 0. 

Thus, by comparison for AP and A/, it is impossible for P — I to have a local minimum at 
Xq. □ 

Corollary 4.3. If xq satisfies -P(xo) = infiP}, the xq € B^. 

Corollary 4.4. If A C M" is the set where P has a second order Taylor Expansion, then 

{xGA: D'^P{x) > 0} C 53(0). 

Now we proceed with the very interesting relationship between inf{P} and inf{F} which is 
imposed by the equation governing F and will be developed as Lemma 14.51 and Proposition 
and 14.71 Lemma 14.51 is an example of the type of information that is available when analyzing 
nonlocal equations which is not available in the local setting. Basically this next lemma says 
that the nonlocal equation prevents, in a uniform fashion, F (and also u) from growing away 
from its minimum faster that lyl'^ in a large set. 

Lemma 4.5. Let < cr < 2 and xq he such that F(xo) = infjRnjF}. Define the set 

A^, := {y : F(xo + y) - F(xo) < /(xq) (4.2) 

and rings with radius r^ = pq2^^ , Rk = Br^{xQ) \ -B^^^^ (where po will be chosen later for the 
proof of Proposition \4. 7^ . If the indices, k* , indicate bad rings for which 

Al.f^Rk' >hRk'\, (4.3) 



then the number of such indices satisfies 



Proof of Lemma \4-5\ Because xq is a location of a global minimum of F which is also a contact 
point with u, we know that M~(F,xo) < /(a^o) by comparison with u. Furthermore, there are 
two things which result from 6T{xq, y) > 0: one, all of the operators appearing in the definition 
of M~ (equation (|1.6|) ) are bounded below by the one corresponding to a matrix {\/n)e®e (for 
some unit vector e), and two, we can neglect as much of the set of integration as we like in the 
evaluation of M~(F,xo). Without loss of generality, we may assume e points in the positive 
direction of the xi-axis. 

Therefore, evaluating (jl.Sp at xq and estimating as suggested above, we have 
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fixo)>M-{T,xo) 
Jm." 

= X [ (2 - a){r{xo + y)- r{xo))yj bl""^'^'" dy + X [ (2 - a)(r(xo - y) - r{xo))yl \y\ 
> A / (2 - c7)(r(xo + y) - r(xo))2/? lyr""""' dy 
>A^ / (2-a)(r(xo + y)-r(xo))y?|y|— '^-^dy 

>AV / (2-a)/(xo)|yry?|ypn-^-2^y^ 



Since 1^4^^ fl-^fe*! ^ l-^fel/^, it can be checked easily that 



/ yflyr'^-^dy > cn Vfc* 



which together with the previous inequahty means that 

/(xo) > Xfixo) 5^(2 - a)c„ 1^^,, fl Rk* 
k* 

> Xfixo) (2 - a)cnOj{n){l - T 

k*+l 

>Xf{xo){2-a)^XCnUJn{#k*) 

Therefore, we conclude 



-n+a-a 
' k*+l 



□ 



Remark 4.6. It is worth comparing Lemma 14.51 to jl8|. Lemma 8.1], where a decay rate of T 
away from its supporting hyperplane of was used instead of the r°" used here. Using r'^ seems 
to give more precise information for which rings the decay estimate can fail. 

Proposition 4.7. Let xq be such that T{xq) = infKnjr} and a € (0,2). Then 



inf{P} = -inf{P} > C(n,A)(-r(xo))^/^ 

R" -B3 



2/(xo)J 



Proof of Proposition \4. 7\ First we note the equality between the infimum of P in all of and 
in i?3 comes from Corollary 14.31 Once we choose an appropriate po from Lemma 14. 5| this 
proposition will follow directly from the definition of P. We would like to restrict the beginning 
radius, po, so that whenever y £ A^g H Rp^, y also satisfies (recall F is always negative) 

T{xo + y) < ^F(xo). 
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This is guaranteed if we let 



which gives for \y\ < po and y € A^^ 



V 2/(xo) ; 



-1. 



r(xo + y)- T{xo) < f{xo) < ^r(xo). 

Also in the estimates to follow, it will be very useful to observe that if {kg} are the collection 
of indices resulting in good rings satisfying 



and N{xq) is the upper bound on the bad indices, {^k*), where the collection {k*} are the ones 
satisfying (I4.3p . then 

„2-(T 



Recall also that 



5^ ^fcg — 

kg 

N{xo) 



k>N{xo) 

2 



(2 - a)XuJn' 



Now we have 



-Pixo) = Ain,2-a) [ -T{xo + 



y) \y\ 



dy 



> A{n,2- a) 

> A{n,2- a) 



1 



r(2;o) \y\ 



-n+2-a 



dy 



> A{n,2- a) 
= A{n,2- a) 
= A{n,2- a) 

> A{n,2- a) 



ln-o))^\Af]R,^ 

kg 

r(xo))C(n)j;: 



kg + 1 



„2-o- 



-r(xo))C(n) Yl d'" 

k>N{xo) 

k>N{xo) 

-r(xo))C(n)/>2-(2-(2-))A^(-o)^(2-(2-^))fc 



fc>0 



-r(xo))C(n)p2— (2-2/{Aa."))(l - 2-(2-'^))-i. 

In the calculations above, the "constant", C{n), was used multiple times as different numbers 
to absorb any extra constants which only depended on the dimension. 

At this point we can collect the various information about pQ, N{xo). Observe that 

lim A(n, 2-a)(l- 2(2-'^))-i = C(n), 

(7— 5>2 

which follows from the fact that A{n, 2 — a) ~ 2 - o- and 2^"'^ rs. l 

as fj — 7> 2. We conclude 

-P{xo) > C(?l)2-2/(A.„)(_r(^^))^(2-<x)_ 
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Recall that po 



/^r(xo)V/" 



and that trivially —P{xo) < — inf P so that 



2/(xo) y 



□ 



-inf P > C(n)2-2/{^^")(-r(2;o)) 

which proves the proposition. 

Remark 4.8. There is a strange but straightforward phenomenon for F resulting from the non- 
local nature of M~. As above, suppose that r(xo) = inf{r}; note xq £ Bi by the assumptions 
on u. By the regularity of F, we can evaluate M~(r,xo) (and also recall F = outside of B^), 
and we see that for an appropriate A with Tr(yl) > A 



f{xo) > M-(r.xo) >{2-a) f 8T{x^,y)y^Ay 



y\ ay 



> {2-a) / 6r{xo,y)y'^Ay\y\ 

J\y\>A 



-n—cr—2 



dy > -r(xo)C(n, A,fj). 



This gives the very strange conclusion that 

/(xo) > -C(n, A,cr)r(xo), 

so long as 0" < 2, which is purely an artifact of the nonlocal nature of the operators (this 
information is not preserved as a — )■ 2). Furthermore, this gives good confirmation that the 
choice of po in Proposition 14.71 will never be too large. 

There is another way of stating Proposition 14.71 which allows for a more intuitive interpreta- 
tion as an interpolation result. 

Corollary 4.9. There exists C{n,X) > such that if cr,T and P are all as above, then: 

||r||oo<c(n)||p||^/2 WfWt''^/' 

Written in this form, the nature of Proposition 14.71 as an interpolation estimate becomes clear. 
To see why such an estimate is to be expected, we prove a closely related result. 

Proposition 4.10 (Interpolation; also as Lemma 4.1 of |43j). Suppose a, F and P are as before, 
then if P is C^'^(IR") we have 

||r||oo<C(n,A)||P||^/2 [Plg-f)/' 

Proof of Proposition \4.lO The proof is tailored after interpolation estimates for Holder norms 
(for similar, but not identical statements, the interested reader should see \28\ Chapter 6, 
Appendix]). 

Since F = (-A)(2-'^)/2p we always have (with A = \A{n, -(2 - a))\ ) 

6P{x,y) 



< A 



|F(x)| = ^ 
5Pix,y) 



M" \y 



n+2-a 



dy 



Bp \y 



n+2-(T 



dy 



+ A 


L 







SP{x,y) 

L,|n+2— (T 



dy 



V/3 > 0. 



(4.4) 



The second integral can be controlled as follows 



6Pix,y) 



Ib^^ \y 



n+2-CT 



dy 



< 2\\P\ 



\y 



n+2-CT 



dy = 2ujn\\P\\L 



ao ran\ 



dr 
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\P\ 



L°°(M") 



2uJn 

2^' 



For the first integral we have 

6P{x,y) 



Bp \y 



n+2-(T 



dy 



< 



\5Pix,y)\ 



sup 



Bp \y\ 



-dy 



(4.5) 



/ \5P{x,y)\ \ a;„ ^ w„ ^ 

= sup —p < [P]cl,l(Rn)—p . 

\\y\<p \y\ J ^ ' ' ^ 

We may now plug inequalities (j4.5p and (j4.6p in (j4.4p to obtain: 



(4.6) 



|r(x)| < Aa;„ 



2-cr' 



|P| 



L°°(R")P 



+ 



O" 



Vp > 0. 



(4.7) 



Then we may pick the p > that minimizes the right hand side of (j4.7p . Since we have a convex 
function of the parameter p we only need to get its critical point: 



dp \2 — a a 



=^ ap 



bp' 



cr-1 



p = \l a/h. 



Finally, putting p = \J ajh gives us the minimum value (basically we used Young's inequality) 

« ^-(2^a)/2^(2^a)/2 , ^^^72^-^/2 ^ ^ ^a/2^(2~<7)/2 

2 - cr a (2 - cr)cr 

Where a = 2||P||2,oo(Kn-) and h = Cn[P]c^A(M.^), going back to (|4.7p we obtain 

irf.ll < 2^+-/^^nCr'^^/^A(n,-(2-a)) (2_.)/2 

The proof is now done, as the term on the right hand side stays uniformly bounded for 
a G (0, 2), which is thanks to the behavior of A{n, —(2 — a)) as o" — > 2. □ 

Before we conclude this section, we must prove one last straightforward feature of the poten- 
tial, P. We must be able to compare the values of P on the boundary of some ball, Br, and 
make sure that the difference between the values on the boundary and the infimum of P are 
still comparable to the infimum itself. 

Lemma 4.11 (Decay of The Potential). There exists a radius, Rcr, depending on a and n such 
that 

- inf {P} + inf {P} >\(- inf {P} 



Proof of Lemma \4-ll Let x G dB^i and xq be such that —P{xq) = — inf(P). We will first 
make a computation with a generic R and then choose the value for R^j at the end. Lemma 
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13.101 tells us that F = in M" \ . Therefore we can estimate the integral 



-P{x) = A{n,2- a 
= A{n,2- a 
= A{n, 2 - a 

= A{n,2- a 

< A{n,2- a 

Pixo) 



/ -T{x + y)\y\ 



-n+2-cT 



dy 



-r{x + y) \y\ dy 



-T {xq + w) \xq — X + w\ dw 



T(xo + w) \w\ 



(\xo-X + w\ 



\w\ 



-n+2-cr 



dw 



X0+WGB3 

n+2-o- 



T(xo + w) \w\ 



-n+2-a 



R 



dw 



6 



Now, choosing so that 



Ra 



-n+2-a 



1 

< - 

- 2 



gives the result. 



□ 



5. det(L'^P) AS AN Integro-Differential Operator on T 



The point of using the potential, P, is that we can convert our fractional order equation into 
a 2nd order equation and appeal to known results. In this way we would like also to recognize 
the operation det(L'^P) as a cr-order operator on T, which one can hope will be comparable to 
M~ by below (proved in Lemma l5.8p . This comparison would then allow us to bring together 
the geometric result involving inf(P) and the envelope property that M~{T) < Cf on the set 
Ku- Fundamental to this pursuit is the formula for the determinant: 



(det(B)) 



l/n 



- inf{Tr(^5) : ^ > and det(A) = 1} whenever B >0. 
n 



The interested reader should consult pH] and for further discussion. The integro-differential 
nature of det(Z)^P) acting on V is developed by making rigorous the formal computation of 
convolving F with the derivatives of the Riesz Kernel. 

Lemma 5.1 (Derivatives of The Potential). Assume that for a.e. x ^Vt CC B^, 



< /icr(F,x) < Ci and W^Wwin) ^ C 



■■1, 



where Ci is allowed to depend on 17. Then P G C^'^{Q,) and the second derivative of P at a.e. 
X (z Q in terms o/F is 

D^P{x) = 



(n + fj - 2){n + fj) 



Ain,2-a)6T{x,y) 



ydSiy 



Id 



\y 



n+a+2 



(n + fj) \y 



dy. 



(5.1) 



Remark 5.2. It is worth pointing out explicitly that Lemma |5 . 1 1 gives enough regularity for our 
purposes later to work with P directly in any estimates involving D^P. 
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Remark 5.3. In particular, our definition of h(j{v,x) gives 

D\v * K,^,) = K{v, x) + ^^^£5^ Id. 

Proof of Lemma \5.1[ We consider as a regularization of P, the convolution of F with a modified 
Riesz Kernel ^^2-0-' where for a > we define 



K2-a{y) if \y\ > a 
Qa{y) if \y\ < oi 



Here Qa{y) is the only quadratic polynomial that makes K2_„ as defined above a C'^'^iW^) 
function across (9-Bq,, 

Q^y) = (-(- + <y-2)\y? + {n + a)a') . 

It will be handy to have the formula for the second derivatives of -fi^2-cr ^-way from dB^ 



A(n,2-cj)|y|-"^'^(n + c7-2) (n + a)2i^-Id if |y| >a 
D'K^.Ay) = { ^ " ^ (5.2) 

_A(n,2-^Kn+^j^ if M < a 

We also use the elementary identities (which follow by reflection and rotation symmetry): 

/ y,yjdS{y) = for i ^ j, / yfdS{y) = - I dS{y) Vr > (5.3) 

JdBr JdBr ^ JdBr 

We denote the regularization of P as Pa := r*K^_^, note that Pa € C^'^(]R") for each a > 0. 
To obtain (j5.ip . we will compute the derivatives using (j5.2p inside of the convolution. Then, 
imitating a classical analytic continuation argument |34t p. 45-47] we will pass the formula for 
D^Pa to the limit a — )• 0"*". First, by the reflection symmetry of -f^2-o- about we may write 



D^Paix) = \ f {T{x + y) + T{x - y)) D^K^^,{y)dy 



(5.4) 



The classical argument relies on computing the average D K2_„ in and out of Ba-, we begin 



by noting that according to (j5.2p 

I D''K^_^{y)dy = -A{n, 2-a){n + a- 2)^Id. 

Next, using integration in spherical coordinates together with ()5.3p it is easy to check that 

f D^K^_,{y)dy = yl(n, 2 - a)(n + a - 2) / jyl""-'^ f (n + - Id) dy 

JBg JBg \ \y\ J 

= A{n, 2 - a)(n + cj - 2) ^ jyl""^'" ( ^'^^j\ d - Id^ dy 



\B I 

A(n, 2 - cj)(n + o- - 2)^Id. 



Thus the integral of D^K2_^ in and out of Ba is the same except for a sign, in other words 

/ D^K^^^{y)dy = Va > 0. 
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This means that the identity (I5.4p stih holds if we add zero to the right hand side as follows 
D^P^{x) = \ f (r(x + y) + T{x - y)) D^Kt,{y)dy - T{x) I D^K^^,{y)dy, 
in other words, 

D^Paix) = \ f ST{x,y)D^Ki_,{y)dy, V a > 0. 
This, together with the pointwise bound 

\D^Ki_^{y)\ < Ain,2- a)C(n)|y|-"-, Vy G W\ a > 0, 
implies for all x € S3 that 

\D'PUx)\<Ain,2-a)Cin) [ ^-^^^^§^dy. (5.5) 

JR" \y\ 

We now make a few remarks to conclude the lemma. First, by construction of K2~<jy 
follows that P°' P uniformly. Second, (j5.5p combined with the assumption that F E H^^^ 
imply that all of the hessians, D^Pa, are uniformly bounded in L'^(Bji) for any Br CC -B3; 
furthermore by dominated convergence, D'^Pa{x) D^P{x) pointwise for all x such that 

mx,y)\ 



\y 



-dy < 00, 



which is a.e. x in B^. Third, after extracting a H^^^ weakly convergent subsequence from P^ 

hoc 



and using uniqueness of limits, we can conclude that P indeed is Hf^^{B-i) and (|5.ip for D^P 



indeed holds. 

Finally, we note that once D^P can be computed a.e. x by (jS.ip . the bound on hcr{r,x) 
immediately translates to a C^'^(O) bound for D^P. 

□ 

Now that we have a formula for D^P,we can compute the rest of the operators needed for 
det{D^P). 

Lemma 5.4 (Integro-Differential Form of Tr{AD^P)). For any positive matrix A and a S (0, 2) 
we define 



A Tr(A) ^ 
A„ — A + ^^I. 



a 



(5.6) 



Then Ty^A^D'^P) defines an elliptic integro- differential operator acting on V . Moreover, it is 
given by the formula 

TiiA^D^P){x) = 



(n + a -2)(n + a) , f , Ay , 

^ JR" \y\ 

Proof of Lemma \5.4\ This is a direct consequence of Lemma 15.11 which says 
Tr{A^D^P){x) = 



(5.7) 



{n + a - 2){n + a) 



A{n,2-a) / 5T{x,y) 



y^A^y 



Tr(A,) 



\y 



n+a+2 



(n + a) \y\ 



dy. 
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Since y'^Acry = y'^ Ay + '^^^ \ y\'^ and -^^^Tv{Aa-) = this integral equals ([57 

□ 

Remark 5.5. Observe that the matrix A^j becomes strictly elliptic, that is, A^ > Aid as soon 
as Tr(^) > cjA. It is thanks to this fact that our estimates apply to equations with kernels 
that are not strictly positive. The authors would like to acknowledge the anonymous referee 
for this important observation as well as their suggestion of introducing the change of variables 
A Af, which significantly helps with the presentation. 

Definition 5.6 (Integro-Differential Form of det(D^P)). The operator V — t- det(D^P)(a;), 
defines a nonlinear integro- differential operator acting on T. At the points, x, where 



r lEi^dy < oo and D^P{x) > 0, 
Jr" \y\ 



this definition agrees with 



V„{r, x) = inf {Ta{T, x): A>0 and det{A) = 1}^ = det{D^P), 
where we use the notation Ta(^,x) := Tr[AD^P{x)). 

Remark 5.7. It is worthwhile to remark on the peculiar form of formulas ()5.ip and ()5.7p . Note 
that the kernel 

y'^By TrjB) 

might take negative values in some directions even for matrices such that i? > 0, so the kernels 
appearing in the formulas are not necessarily positive. The kernels will be positive for matrices 
of the form A^^ := A + '^'"^"^^ J, ^ > which is why they were introduced in Lemma 15. 4[ 
However, matrices which are not of this form appear in the formula for the operator V^, which 
means that T)a-{T,x) = det{D^P(x)) is not an elliptic operator for a general F (not even if 
Ea(X) > 0). Despite these issues, serves a crucial purpose in the proof of Theorem 11.31 as 
suggested by the the next two Lemmas. 



f \6T(x,y)\ 

Lemma 5.8 (Nonlocal Determinant - Minimal Operator Ordering). If / ^^^^ — dy < oo 

and Efj{T,x) > 0, then whenever D'^P{x) > 0, it also holds that 



V„{T,x)<^{M-{T,x)r. 



Proof of Lemma \5.8[ Thanks to the assumption of D'^P{x) > 0, we have 

(P.(F,x))^/"< inf |-Xi(F,a 
dct(A)=i [ n 



< inf I-TaA'^^^ 
det{A^)=l,A>0 [n 



AId<£A>0 i n det(A,) i/n ^"^^ 
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Now, by definition (see (15.60 ) of A^j in we have the set inclusion (see also Remark 15. 5p 

{A^ : A> 0, Aa > Aid, A} D {A^ : A> 0, Tr(^) > aX} 
Which together with the previous inequalities shows that 

{VAr,x)f- < ^ inf I^A{n,2-a)[ 6T{x,y)-y^dy] 

A>0,Tr{A)>aX [ A 7]gn \y\ J 

and the Lemma is proved. □ 

Remark 5.9. It is worth remarking that Lemma 15.81 is one of the three key reasons for a good 
definition of the envelope, T. In fact, this ordering between Vu and M~ is not to be expected for 
a general function, v, which does not necessarily satisfy Efj(v) > 0- just as in the second order 
setting between the Monge-Ampere and Minimal Pucci operators within the class of convex 
functions. 

The restriction to points x where D'^P{x) > is not too severe, as later on we will reduce to 
those X for which P coincides with its convex envelope, and there we shall always have D^P > 0. 
Very much related to this fact is the following crucial Lemma which will also guarantee that 
we will only need to test the equation at points where u = T. 

Lemma 5.10 ({det(-D^P) > 0} is contained in the contact set). If for a.e. x, 

f \ST{x,y)\ ^ 
JM" \y\ 

and r 7^ 0, then for a.e. xq such that r(xo) 7^ u{xq), we can find some direction r such that 
Ptt{xq) < 0. In particular, 



({dei{D^P) > 0} [^{D^P > 0}) C{u = r}. 



Proof of Lemma \5.1(K Lemma 13.111 plus standard techniques in viscosity solutions tell us that 
at almost every such xq, EaiT^x^) = 0. Hence for each of these xq there exists a direction, r, 
such that 

{ha)ij{T,Xo)TiTj = 

Hence we can plug r into the formula (15. ip for Ptt and get (recall TjTj = |tP = 1) 

D tu ir w n (' '^r(xo,y) [ 8T{xQ,y) 

Ptt = [h„{r,xo))ijTiTj - Cn,a / , dy = -Cn,a / , dy < 0. 

JR" \y\ JR" \y\ 

In this last inequality we have used the fact that T solves Ea(r) > everywhere, and so because 
— (— A)'^/'^r = Tr(/io-(r)), we have the correct sign for the leftover term. This concludes the 
first assertion of the lemma. 

To finish with the set inclusion, we note from the CoroUarv 14.41 that {D^P > 0} C -B3, and 
the fact that in B3 C\{u 7^ F}, F solves £'o-(F, x) = (see (jl2.4p ). The set inclusion then follows 
by the first assertion of the lemma. □ 



ABP Estimates For Integro-Difforcntial Equations 



25 



6. Proof of The Main Theorem 

We are now in a position to prove the main theorem, Theorem 11.31 The main step left to 
prove is an estimate on inf(P) which involves a measure theoretic norm of /. This will come 
from the geometric set-up of the second order ABP estimate for P, thanks to the various results 
of section [5l 

Proof of Theorem \1.3[ The theorem will first be proved under the assumption that u is C^'^ 
from above, and then we will remove this restriction at the end of the proof using the standard 
method of inf-convolution of u. Therefore, assume u is C^'^ from above. 

As mentioned in the lead-up to this proof, we will be applying the Aleksandrov Estimate 
for convex functions to the convex envelope of P, but a crucial factor to be determined is the 
domain in which the argument will be applied. We must work with an appropriately large 
radius according to Lemma 14.111 in order to make sure that the difference of the values of P 
on dB^^ and inf(P) are still comparable to infiRn(P). Therefore we are working in Bp>^ from 
Lemma 14.111 

Let P^^ be the convex envelope of P in Br^. It is essential to note a very important 
consequence of Lemma 14.21 is that the contact set, {x : P{x) = P^^{x)}, is contained within 
B3. Furthermore, since P £ Cj^liBa), we know that P"^^ G CI;1{BrJ (see [lH Chapter 3]). 
The Aleksandrov Estimate for convex functions (see [21^ Chapter 3]) then implies 

- inf {P} < (- inf {P}] +C{n) ( [ det{D^P^^{x))dx\ 
< (-^mf {P}^ +C{n) ^ det{D^P{x))dx^ 



(- inf {P}] +C{n) ( / V^{r,x)dx 

V 9^R<T J \j!P=pCE} 



1/n 



We note that the constant C(n) in fact depends upon R^j, but it is uniformly bounded from 
below and above for our R^^, and so we simply denote it as a dimensional constant. 

Thanks to the comparison between D„ and M~ , Lemmas 15.81 and 15.101 and the comparison 
principle between F and u for M~ on all combine to give 



mf{P} < (-^mf {P}) +C(n) P.(F,x)dx) 



l/n 

1/n 
l/n 



We note that we have also used the C^'^ from above nature of u to evaluate (jl.Sp classically. 
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In order to conclude the theorem in the case that u is C^'^ from above we must now recall 
how inf(P) relates to inf(r). We can now combine the results of Proposition 14.71 and Lemma 
14.111 with (j6.ip to obtain (with obvious abuse of the dimensional constant, C{n)): 

c,„,.)(-r,.o),^/"(^) ^ (4 r..^"" 

which after rearranging gives 

-inf{r} < ^(||/||L^(^.))(^-'^)/^(||/lk.(K.)r/^. 

Now to finish the theorem for a generic u, we work with the inf-convolution of u, 

Ue{x) = inf {u{y) + - \x - yf} 

k-s/|<v^ll'"lloo £ 

(see [m Appendix A], |25l Appendix], [301 Equations (14), (15)], or [35] for definitions and 
[13\ Appendix A] for the analogous argument regarding the ABP in the second order setting) 
The main properties we will use are that is C^'^ from above in B^, that is increasing to 
u, and that solves the equation (jl.Sp with / replaced by (see [13\ Lemma A. 3] and [181 
Proposition 5.5]) given as 

fix) = sup {/(y)}. 

For Fg, P^, and K^^ denoting all the corresponding operations using in place of u, ()6.ip 
becomes 

C(n) ' " ^ ^ 




So long as limsup^ /^(x) < f{x), which is given by continuity of / in this case, we can conclude 
the estimate by letting e — t- and using Lemma 13.161 □ 

Remark 6.1. For future reference, we would like to collect a very important fact which is a 
cornerstone of the proof of Theorem 11.31 It is the relationship between inf(F) and the integral 
of V„(T), 

-inf(r)<^f/ V^{r,x)dx) (6.2) 



A 

which holds whenever 

-dy < oo 

Jr" \y\" ' " 

for a.e. x. 



mx,y) \ 



7. Fractional Monge-Ampere Type Operators 

The Monge-Ampere type operator T)„{r,x) is crucial to our proof of Theorem 11.31 ^-s it 
allows us to borrow the divergence structure of the standard Monge-Ampere which is essential 
for the classical ABP Theorem. One may take the point of view that T„(x) = VP^ is a nonlocal 
gradient map given by u, and that is just the Jacobian of the map T^. Then, part of the 
effort in our proof (Section U]) has been relating the infimum of F to the size of the image of the 
map Ty-{x). 
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We insist in referring to T>^ as a "Monge- Ampere type operator" and not just the "Monge- 
Ampere operator", as it is not clear that would be a definite, canonical analogue of this 
operator for integro-differential equations. Note also that as Pu is not necessarily convex the 
map VPu might not be monotone, and in fact the operator T>fj[u) might fail to be elliptic even 
on the class of functions satisfying > (for the potential P might not be convex). Another, 
perhaps more natural candidate for an "integro-differential Monge-Ampere" is given by the 
operator 



Where n„ may be chosen to he n^j = n or at least such that n„ ^ n as a ^ 2. Unlike 
Per, this second operator is affine invariant and it can be checked easily that it is well defined 
and (degenerate) elliptic in the class of functions which are subsolutions of = 0. The main 
disadvantage of P* with respect to 2?o- is that it is hard to relate directly the size D*r to the 
infimum of F. This relationship could be quantified if P* had a connection to some sort of 
gradient map, which is (informally speaking) how the divergence structure of Monge-Ampere 
contributes to the classical ABP (recall V^j = det(DT„)). 

The operator as defined in (|7.ip would be a natural operator to consider in the class of 
equations given by 



A posteriori, one can see the reasons for the proof of Theorem [L3] being restricted to a much 
smaller class of equations than those appearing in |18] . There must be a balance to ensure 
each one of Lemmas 15.81 15.101 and (16. 2p hold. In principle, one could create a new notion of 
determinant altogether, such that those three meta-lemmas would still be valid for much richer 
families of kernels. The difficulty comes from the fact that considering a larger family of kernels 
makes the operator M~ more extremal, and thus the question of finding a "geometric" operator 
that is comparable M~ becomes much harder to tackle. In conclusion, all these issues underline 
the lack of geometric equations for integro-differential operators. 



In the series of works, [16], [T7j, and [18], all of the results were obtained in a fashion in 
which they are preserved as cr —t- 2 and recover the corresponding results already proved for 
second order equations. This has led to a unified picture of the second order and the fractional 
order (i.e. nonlocal) theories. We adopt this view in the current work, and take this section to 
discuss explicitly how our result relates to the relevant second order theory. 

8.1. Recovery of A Second Order Envelope as cr — )• 2. One may expect that as o" ^ 2 the 
envelopes should behave more and more like the convex envelope of u. This is almost the 
case at least whenever u is C^'^ from above, the discrepancy arises only because of the behavior 
of E'er as fj —T- 2 (see Section [3]) . 

Proposition 8.1. Assume u is C^'^ from above, then as a ^ 2 the envelopes To- converge 
uniformly to a function T2 which solves the obstacle problem 




(7.1) 




where A^ = A, AI < < AI 



(7.2) 



8. Important Limits As a ^ 2 



m.m{u{x) -r2{x),L{r2,x)} = y x £ B3, r2(x) = Vx G dB3. 
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where L{v, x) is the second order fully non-linear elliptic operator defined by 

L{v,x) := \i{D'^v{x)) + \i:^v{x). 
In particular, Fq lies in between u and the convex envelope of u in B3. 

Proof of Proposition Here we will use the operator, E*, which was introduced in Remark 
13.81 Let To- and F* be respectively the envelopes of u made using Definition l3.6l with respectively 
the operators and E*. We note here that because E*^^ > in i?3, then it also holds that 
-^'o-T* > in i?3 as well. Moreover, T* satisfies T* < u by construction. Hence because To- is 
the supremum of such functions, we see that 

r* < r 

Due to the C^'^ assumption on u we can use Lemma 13.111 to conclude that both families, 
{ro-}cr and {r*}o- are uniformly C^'^ from above. In particular there exists a constant C > 
independent of a such that 

_(_A)'^/2r^ < C and - (-A)'"/2p* < 

Since we also have by construction both 

_(_A)'^/2r^ >0 and -{-Ay/^Tl>0, 

we conclude that both \\{-AY/^T^\\lc^^Bs) < C and \\{-Ay/^Tl\\Loo^Bs) < C. Since each 
and r* are identically zero outside B3 we may use (for instance) the Poisson kernel for (— A)'^ 
in [31] to conclude both T„ and L* are Holder continuous in M" and uniformly in a. 

Therefore, from each sequence Ufc ^ 2 we may select a subsequence that converges uniformly 
in M" to some functions r2 and Fg. 

We note, for any smooth v with enough decay at infinity, the convergence of E„{v,x) as 
o" ^ 2 to Xi{D'^v{x)) + ^Av{x) (see Remark 15.31 and Proposition 18. 2p and the convergence 
of E*{v,x) to Xi{D'^v{x)). Thus, by the stability of viscosity solutions under uniform limits, 
we conclude r2 and r2 are the unique viscosity solutions to their respective limiting obstacle 
problems, and hence all subsequential limits converge to the same function. But in the limit 
r2 is the convex envelope of u, which concludes the proposition. □ 

We proved the proposition assuming u is C^'^ from above for the sake of simplicity, but as it 
is reasonable to expect the envelopes Tg to be somewhat regular regardless of u, it is likely this 
convergence holds under much more general circumstances. 

8.2. Derivatives of P Converge To Derivatives of v. We now revisit the peculiar form 
of (15. ip and (15. TP mentioned in Remark 15.71 We give further justification of formula (15. ID by 
checking with a direct computation that when P := (— A)"^^""")/^?; then as u — >■ 2~ it gives 
back the Hessian of v, assuming v is smooth enough. This is not surprising since P, being 
defined as (— A)~(^~'^-*/^u converges to v uniformly as cr ^ 2~, so if f is smooth enough we can 
already conclude D^P — )• D'^v, in this regard the proposition below is only complementary. 

Proposition 8.2. If v is a fixed function which is C^'^ from above in B3 and P is its Riesz 
potential j^^. j[ ), then 

lim D'^P{x) = D'^v{x), V X € S3. 

0--5.2- 
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Proof of Proposition \8.2l The proof has two broad steps. 

Step 1) We shah compute a formula for Hm Px^x ■ To do this, we estimate the integral on 



(T-s>2- 



the right hand side of ()5.ip by breaking it in two parts. Fix e > 0, observe that 

ViVj Sij 



{2 -a) 6vix,y) 
B? 



\y 



\n+(T+2 



{n + a) \y\ 



dy 



I 1 /■+00 ^ ; 



dt 



The right hand side vanishes as <t — >■ 2 (for e > fixed), thus 



lim Pr 



n(n + 2) 



XIX] 



lim 



A{n,2- (7 



^ ^v{x,y) 

|y|n+(T 



ViVj 



Si 



{n + a) 



dy 



.11 



2 cr^2-7B,(o) \yy'" i\y\ 

This holds for all e > (which makes clear we are getting a local operator). To estimate the 
integral inside -B^ recall that v is C^'^(x), which means that Vv{x) and D^v{x) exist in the 
sense that as \y\ — > 

v{x + y) = v{x) + Vv{x) ■ y + \y'^ {D'^v{x))y + o(|yp) 



In particular, setting y = y/ \y\ 

Sv{x,y) 

Using this expansion in (18. ip we get furthei0 

n(n + 2 



f iD'vix))y + o{l) , |y| ^0 



lim P-r 



lim <(2-a) 

CT^.2 



a^2- ' 2a;, 
We may also use polar coordinates to see that 



y'^ {D'^v{x))y 



B, \y 



n-2+o- 



ViVj 



n + a 



dy 



y'^{D'^v{x))y 



Be \y 



n-2+(T 



ViVj 



n + a 



dy 



e^{D'^v{x))e 



1 



,2-0- 



2-0- 

This gives us the formula 

n{n + 2) 



e^{D'^v{x))e 



n + a 



n + a 
dS{e) 



j.n—1 



Z2+^dS{e)dt 



lim Pr 



0-^2- 



2a;„, 



eiCj - ^^-^-^Jjj ) e^{D'^v{x))e dS{e) 



i.2) 



Step 2) All there is left to show is that the expression on the right in (18. 2p always gives back 
the ij entry of D'^v{x). By rotation invariance we may assume without loss of generality that 

n 

D'^v{x) is diagonal, in which case we have {D'^v{x))e = ^ vii{x)ef for all e, in other words: 

1=1 



n(n + 2) 

lim_ Px^xj = — 5 tntgn-i [ eiCj 



n + 2 



ef dS{e)vxixi{x) 



^&\so recall that A{n,2 — a) / {2 — a) u)^^ as ct — > 2 
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Note that by reflection symmetry the integral of eiCjef must be zero for all I when i ^ j, given 
that D'^v{x) is diagonal this means lim Pxix = = Vx^x ix) for i ^ j. We are left to consider 

the case i = j, to fix ideas we do it for i = 1. Then 

lim Pxixi = — — — — ^ f I e? I ef dS(e)vx,x,(x) 

Due to rotation invariance this can be rewritten as 

^ [/,..-. - M ''"<'>°"'<^> + (/..,-. ('-^^ - ^) g"-'<^' 

One may compute explicitly the integrals above and get 
/ e^d5(y) = ^^, / e\e]dS{y)- 




Therefore 

4 - ^1 dS{e) = [ ( elel - ) dS{e) = 

^ n + 2j ^ ' n(n + 2) y^n-i V n + 2 ' ^ ^ 

Plugging this in the last expression for lim Px-^^x^i^) obtain 



lim PxrxA^) = 7, ( TTi^'^^i^i + I - ux^xi 



n{n + 2) / 2ijJn 

This proves the proposition. □ 



Remark 8.3. The spherical integrals above are not hard to compute, they follow from counting 
the different terms appearing in the trivial identities: 



/ el + ... + eldS{e)= I {ej + ... + elf dSie) = [ dSit 
Js"-'^ JS"-i ./5"-i 



UJr, 



and from the relation /^„_i ef e2(i5'(e) = | Js^-i e.fdS{e) which is a standard computation for 
n = 2 and can be pushed for all n via induction (integrating along slices of the sphere and 
rescaling the lower dimensional formula in the inductive step). 

8.3. Operators covered in the limit cj — > 2^. We now characterize the second order opera- 
tors that can be obtained as a limit of fractional order operators of the form (jl.3p . It is rather 
surprising that many elliptic linear operators of second order cannot be obtained as limits of 
operators of the form (jl.Sp . as opposed to what was expected (see discussion in [181 Section 3]). 

Proposition 8.4. Given a symmetric n x n matrix M and v € Cf^^{W^) n L°°(]R") we define 

LmAv, x) = TriM ■ K{v, x)) = (r^ + ^ - 2)(n + a) ^ _ f y^V^^^y^ 
Then, for any x € M" we have 

lim LmA^^x) = Tr ({M - • D\{x) 
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Moreover, an operator of the form L{v, x) can be obtained as a limit of the 

operators L M,a for some AI > if and only if 

A - ^Id > 0. (8.3) 

Proof. This proposition is merely and application of Lemma |5.H which in our case says that, 
if we consider the function = (—A) ^ 2 >v then 

\y\ ^ ' 

\y\ 



In particular, we get 

-IdMpH) 2 - .) / "^dy = -1AP(.). 

jr" \y\ 

Then, again by Lemma |5. II and the definition of h^{v,x), we have 

h^{v,x) = D^P{x) + ^^U. 

Thus, 

LmAv,^) = Tr (m . (D^P^^ix) + IlliZ^Id)' 



= Tt{M ■ D^P^{x)) + ^Tt{M)Tr{D'^P^{x)) 
= Tr((A/ + 2^id).Z)2p-(^)^ 

Since v £ Ci^^{W^) we have D^P^ — D'^v{x) for every x as cr — )■ 2~, and 

lim_ LmA'"^ ^) = Tr(A • D^v{x)), A := M + 2^Id. 

Now, we invert the relation defining A in terms of M, it is worth noticing that this relation is 
nothing else but A = M„ when o" = 2, as defined previously in (j5.6p . Taking the trace on both 
sides gives us Ti{A) = (1 + §)Tr(M), so that M = A - ^^Id = A - ^^Id, from where it 
follows that M > if and only if ^ — -^^Id > 0, and the proposition is proved. □ 

Remark 8.5. Note that for A = Aid (A > 0) we have A — -^q^Id = A^^Id > 0. In particular, 
in the space of matrices there is a tubular neighborhood of this ray on which (j8.3p holds. 



Remark 8.6. As it was natural to expect, the operator L{v,x) appearing in Proposition 18.21 
(once linearized about some v) satisfies (|8.3p . 

9. Comparison Theorems Related to Theorem 11.31 

In this section we collect various results that are either direct applications of Theorem 11.31 
or straightforward modifications of its proof. Each of the results here are stated without proof, 
and we simply mention some of the modifications. First we mention in Theorem 19 . 1 1 the analog 
of Theorem 11.31 to more general domains than just Bi, and in Theorem 19.21 the applications to 
the comparison principle for (jl.Sp . 
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Theorem 9.1. Let equation 111.8]) be set in a general bounded, connected, domain, D instead 
of in Bi . Then 

-infM < ^diam(i?)(||/||i.(^„))(2--)/2(||/||^„(^^))'^/2. 

We note that to modify the proof of Theorem [T3] to incorporate -D, one simply needs to modify 
the domain of truncation in the definition of r„ , Definition 13.61 ^-^id also the selection of the 
radius, R^, from Lemma l4.11i In particular, if we define the set D3 := {x : distance{x, D) < 
3}, then we are concerned with the obstacle problem in D3, using uId in Definition 13. 6i 
Furthermore, the ball -63 is no longer used in the proof of Lemma 14.111 but instead is 
chosen large enough to compensate for the size of instead of the size of B3. 

An immediate consequence of the fact that the difference between a subsolution and a su- 
persolution with solve the minimal equation for M~ is that we get a comparison theorem for 
subsolutions and supersolutions of equations with the same operator, F, but different right 
hand sides. 

Theorem 9.2. Suppose that F is in the elliptic family for M~ and M^. Letu andv be bounded 
and respectively use subsolution and Isc supersolution of 

{F{u,x) > f{x) in D 
u = uo on M" \ D 

and 

{F{v,x) < g{x) in D 
v = vo on \ D. 

then 

sup{n — v} 

D 

< sup {uo - vo} + ^diam(Z?)(||(/ - g)-\\L^iK._^))^'~''^/\\\{f - grUnK...))"/' ■ 

10. Special Cases of the Regularity Theory of Caffarelli and Silvestre 

In this section we show how Theorem 1 1 . 3 1 can be used to prove the standard estimates for 
viscosity solutions of (|1.8p (at least for a special family of operators). None of these results are 
new, as such estimates have already been proved in \18\ Section 10] for very general nonlinear 
integro-differential equations. The estimate constitutes the backbone of the regularity theory 
for fully nonlinear elliptic equations; it is the key fact behind the Harnack inequality of Krylov- 
Safonov, the Evans-Krylov theorem and the respective Caffarelli-Silvestre theorems for nonlocal 
equations (again, see |21j and |18|). Our purpose in revisiting this part of the theory in our 
case is basically illustrating the uses of Theorem 11.31 in particular Theorem 11.31 allows us to 
do essentially the same proof of the L*^ bound for second order equations used in |2H Lemma 
4.6] (compare to the different method needed in jl8|. Section 10]). Furthermore, the estimates 
obtained are uniform in the order of the equation, recovering the second order theory as o" — )> 2, 
which was already done in jlSj 

We state without proof the following proposition regarding the construction of a special 
barrier function. It follows by an argument similar to that used in Lemma 9.1 of |18) . 
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Proposition 10.1. Given < A < A and ctq € (0,2) there exist constants Co,M > and a 
C^'^ function r]{x) : M" — )• M such that 

(1) supp r] C -62^(0) 

(2) T] < —2 in Q3 and ||?7||oo ^ M 

(3) For every a > ao we have M^{ri,x) < Cq^ everywhere where ^ is a continuous function 
with support inside -61/4(0) and suc/i that < ^ < 1. 

With this special function in hand, Theorem 11.31 allows (just as in the second order theory) 
to control the average size of a supersolution in terms of its value at a point. This point to 
average estimate implies a weak L'^ estimate for supersolutions which is the key step in the 
proof of the Krylov-Safonov and Caffarelli-Silvestre regularity theorems. 

Lemma 10.2. Given n, A, A and ctq such that < A < A and (Tq € (0,2), one can find positive 
constants M > < 1 and 5q such that if u satisfies: 

(1) li > in W 

(2) miu < 1 

Qa 

(3) M-{u,x) < f in (for some a > ao), and Wfh^iQ^^J < So, < 1- 
Then we have the bound 

\{x e Qi : u{x) < M}\ > fi\Qi\. (10.1) 

Proof of Lemma \10.2l Consider the function w = u + r] {1] as in Proposition 110. ip , then w 
satisfies (in the viscosity sense) 

(M-{w,x)<f{x) + CoC ini?4^ 
\w>0 on M"\S4^. 

Moreover — inf {w} > 1 since u < 1 somewhere in and r] < —2 everywhere in Q-^. In this 
Q3 

situation. Theorem 11.31 (rescaled to the ball B^^) tell us that 

1 < CA-" (1 + Coi^-"^/^ (11/ + coiu^^K^^r'^ 

where we recall K^^ = {x € B^^ : w{x) = r^(2;)}, then 

C"2/'^A2«/-(i + Co)"(2-)/- < 11/ + CoC||i„(^„). 

Hence 

fj~2/ay2n/a^-^ + Cq)"'^'^^"^/'' < 5q + Cq\K^ R ^1/4!^/". (10.3) 

One then sees that picking 5q universally small one gets for a universal /x € (0, 1) the lower 
bound 

fi\Qi\<\K^r\Bi/i\<\K^nQi\. (10.4) 

Now X G implies in particular that < therefore u < —t] < M, where M = sup \r]\. Since 
Kw D Qi C {x Qi : u < M} this last inequality proves the lemma. □ 

As done in |2H Chapter 4, Section 4.2] one can use the Calderon-Zygmund decomposition 
and 110.21 to prove the weak L*^ estimates we previously mentioned. 
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Theorem 10.3 (weak L*^ estimate). Let u > in M" be a supersolution of M {u,x) < f in 
Bi (for a > (Tq) and such that ti(0) < 1. Suppose that 

||/||oo<l, ||/IUn(B,) < 5o, fTe((Jo,2). 

Then there are universal constants C, Jq oi'^d e (i.e. determined by n, uo, A and A) such that for 
all t > we have 

\{u>t}nBi/2\<Ct-^. (10.5) 

With this lemma in hand, that viscosity solutions to (jl.Sp are Holder continuous or that they 
satisfy a Harnack inequality follows by standard arguments that can be found in |2H Theorem 
4.3 and Proposition 4.10] for second order equations and in |18| Theorem 12.1 and Theorem 
11.1] for integro-differential equations. 

11. Applications and Open Problems 

Here we would like to make a few remarks about further research directions and open ques- 
tions where we anticipate Theorem 11.31 could be useful, some of which are in obvious analogy 
to the second order theory. Accordingly, the discussion below is only suggestive but we include 
it with the hope of stimulating further work. 

Stochastic homogenization. As pointed out in the introduction. Theorem 11.31 will play an 
important role in the homogenization of stationary ergodic families of equations within the 
ellipticity class governed by M~ , and this will be presented in |41j . 

up^cr,£5? ^j^jj up^CT,p5? estimates. These would be analogous to the W^'"^ estimates of Caffarelli 
[211 Chapter 7] and Lin [37] and subsequently the W'^'P theory of Caffarelli |2D]. Such estimates 
would allow use of the regularity theory of [HI [T71 [16] to obtain W^'^ regularity for viscosity 
solutions of ()1.5p in terms of the norm of the right hand side /. It is important to remark 
that such estimates are not yet available for any kind of fully nonlinear equation of fractional 
order. 



More general equations. Comparison and regularity theory for more general equations where 
their ellipticity is considered on choices of (with varying scopes of generality) such as: 



LA{v,x) = {2-a) 6u{x,y)a{x,y)\yr''"'dy (11.1) 
where a(x, y) is homogeneous of degree in y, 

La{v, x) = {2-a) [ 5u{x, y){y^ A{x)y)n{dy) (11.2) 



where n{dy) is comparable to \y\ ^ ^ dy only in some subset of M" (indicated to the authors 

in m), 

LA{v,x) = {2-a) [ 5u{x,y)^^^dy, X<a<A (11.3) 
JR" \y\ 

which corresponds to the family in [18], and 

La{v,x) = {2- a) 6u{x,y)n{x,dy) (11.4) 
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where n(x,dy) (for each x) is again a measure comparable to \y\ dy only in some subsets 
of M" (for example, along certain directions through the origin). 

The "right" nonlocal Monge- Ampere equation. It is worthwhile to find out whether 
another notion of nonlocal determinant can be found which is both extremal and carries a geo- 
metric interpretation (or "divergence structure") such that it achieves the key features listed 
in Section I2.5l - (j2.8p . ()2.9p . and (j2.10p - for a more general family of equations, such as those 
considered in [TS] . 



12. Appendix- Concave Nonlocal Obstacle Problems 

In this section, we develop some regularity results for solutions to obstacle problems with 
concave nonlocal operators. We believe these results may be of independent interest and have 
hence chosen to present them in a generic form, but the example we actually use in this work 
is for the operator, E^. The results in this section are nonlocal analogs to the results of |12j 
and [36]; the interested reader can also find them presented in the book [331 Chapter IV]. 

First we list some notation and assumptions for this section: 

(1) V : -B3 ^ M is the obstacle function with ■0 < G ^3 and V > G \ S3. 

(2) F is a degenerate elliptic, translation invariant, nonlocal operator which is concave 
with respect to v in the sense that there is a collection of linear cr-order operators, C, 
and F can be represented as 

F(v,x) = mi{L(v,x)}. (12.1) 
Lec 

(3) F satisfies a non-degeneracy type condition as 

- {-A)''/^v{x) > F{v,x). (12.2) 

(4) j3^{s) = /?(-) where /? : R ^ M is convex, non-negative, decreasing, /3(0) > sup^g^^ -^(^^ x) 
and (3{s) = for s > 1. 

(5) CC ^3- 

The main function of interest for this section is the solution to the obstacle problem, 
solving 

V = sup{w : F{w,x) > in S3 and w <i^lB; in M"}. (12.3) 
It is standard that v is the unique viscosity solution of the equation 



(12.4) 



{min.{F{v, x),u — v} = ^ in -B3 
v = Q on M" \ S3. 

One of the main tools for investigating ()12.3p is the penalized equation, which also admits a 
unique viscosity solution: 



F{v',x)=l3e{i:-v') in S3 

= ^ on M" \ S3. ^ ' ' 
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Finally, in order to make some arguments easier by working with classical solutions, we will use 
the regularized, penalized equation: 

(F{v'P,x) + pAv'P = f5,{i;-v'f) inBs 

[v^f = on M"\53. ^ ' ^ 

The main result of this section says that the solutions of these obstacle problems not only 
inherit the regularity from above of the obstacle, but in fact they are even regularizing in the 
sense that they obey interior regularity from below as well, which may in fact be more regular 
than the obstacle, ^. 

Proposition 12.1. If^l' satisfies classically —{—A)°'^'^'tp(-) < C and Q CC -B3 then there exists 
a constant, Ci, depending only on C, a, n, dist{Q,dB-j), and such that for a.e. x G 0, 

V satisfies 

\y\ 

and 

< -{-Ay/'^v{x) < Ci and ||r||^a(f^) < Ci. 

Remark 12.2. It is worth pointing out that the same results of this section hold mutatis mu- 
tandi if one replaces the operator — (— A)*^/^ by any other linear, translation invariant, a-order 
operator (and its corresponding kernel, K) which has similar regularizing properties to those 
of — (— A)'^/^ used in the proof of Lemma ll2.51 

In fact, if i{j has enough regularity such that for all of the translation invariant L in the ellipticity 
class of [m Section 3] (the operators governed by the extremal operator of (|2.4p ) it holds that 
x) < C, then a much stronger result holds. The proof is a simple calculation, and is found 
in Section 7], but we include it here for completeness. 

Corollary 12.3. Let L he the collection of linear operators, L, of the form 

L{w,x) = [ 8w{x,y)K{y)dy where K is measurable and -j— < K{y) < -j—. 

Jr" \y\ \y\ 

If ip satisfies L^ipjx) < C for all L G C, then there exists a constant, Ci, which depends only 
on C, a, n, \, A, dist{0,,dB3), and 11-011/, 00 such that for a.e. x (z 0, CC B^, 



\y\ 



\y\ 

Proof of Corollary \1 2. 3[ Proposition II 2.11 (modified for L), says that for all L € 

< L{v,x) < Ci, 
which in particular implies (see [181 Section 3]) 

A {6v{x, y))'^ - X {6v{x, y))' 



\y\ 



-dy < Ci. 



On the other hand by ()12.2p and (jl2.3p . we have that 

0<-(-Ar/\.(.)=/ '^,y, 

JK" \y\ 
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and hence multiplying through by —A 

-X{6v{x,y))^ + X{6v{x,y)) 



dy <0. 

y\ 

Appropriately adding together these estimates gives and estimate on {6v{x,y))^ \y\~^~'^ i which 
can be plugged back into the estimate < — (— A)°'/^f (x) to conclude the Corollary. □ 

The proof of Proposition 112.11 is broken up into multiple pieces which are worth recording in 
their own right as Lemmas 112.41 and 112.51 Once these results are proved, Proposition follows 
immediately. 

Lemma 12.4. Let s > be fixed with s < dist{Q, B^), and let the truncated kernel and operator 
he defined as 

K{y) ■.= \y\~''~'' l\y\<s and L{u,x) := / 5u{x,y)K{y)dy. (12.7) 



Let v^P he the unique smooth solution to the regularized and penalized Dirichlet prohlem \12. 6]) 
and the unique solution of \12. 5j) . Then v^^ — )■ f ^ as p ^ locally uniformly, v^^ satisfies 
the estimate 

ma^{L{v^P,x)} < max{L(V', x)}, (12.8) 
and v'^ satisfies in the viscosity sense 

L{v^,x) < max{L(?/',x)}. (12.9) 

Proof of Lemma 12. 4\ The whole point of working with v^^ is simply to have a classical solution 
to manipulate more easily. By the concavity of -F, as (jl2.ip . and the linearity of L, the functions 
w := L{v^P) satisfy in the equation 

F{w,x) + /\w{x) > L{/3s{i; -v'P),x). (12.10) 

The main point of choosing /3 to be convex is the fact that this implies 

L(/3,(^-0,^) > -m\L^L{{i;-v'P),x). (12.11) 

To see this, we note that convexity gives 

f3, ii^P - v'P){x + y)) + I3e ii^P - v'P)ix - y)) > 2/3, f ^(^ - v'P)ix + y) + i(^ - v'P){x - y) 



X 2 
which after subtracting 2{ip — v^p){x) and estimating with the Lipschitz norm of gives 

6f3, o (V; - v'P){x, y) > -\%\\oo m " v''){x. v)) ■ 
Thus (jl2.1ip follows. Therefore, plugging back into (jl2.10p gives 

F{w,x) + /^w{x) > -\\^,\\l^L{{'^-v'P),x), 
which can be evaluated at Xmax such that t(;(xinax) = max{w} to give 

> F('U;,Xmax) + Au;(Xmax) > " || /^e || L°° ^((V' " ^'^'') , 2;max) ■ 

We note the use of the simple fact that at a maximum point, both F{w) and Aw are nonpositive. 
Therefore, multiplying through by — ||/3e||L°° gives (jl2.8p . 

The fact that (|12.5p admits a unique viscosity solution, combined with the stability of vis- 
cosity solutions of (|12.8p immediately gives via standard techniques that v^p as p ^ 
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locally uniformly (see any combination of [5], [181 Sections 4,5], [25] )■ Finally, the stability 
of the inequality (|12.8p with respect to local uniform limits immediately gives (|12.9p . and we 
conclude the proof of the lemma. □ 

Lemma 12.5. Let v be the unique solution of ^12.4\ )- There is a subsequence, still denoted as 
, such that ^ v locally uniformly in B3 as well as weakly in H'^ as e —?■ and 

< -i-Ay/^v'{x) < Ci, 

in the viscosity sense and where Ci depends only on max^gl— (— A)'^/^V'}! IIV'II-L°°; ^■ 

Proof of Lemma \12.5[ We start with the observation that constants are subsolutions of F = 0, 
and so V > inf^gl'i/'}. Furthermore, by the definition of L and (jl2.9p . we have 

_(_A)-/V=/ 5v'{x,y)K{y)dy+ f 6v' {x,y)K{y)dy 

J\y\<s J\y\>s 



< max 

< max 

Bs 



[L{^)]+C{s,n,a)\mL^ 
|_(_A)-/2V;}+C(s,n,a)||V'||L-, 



where C{s,n,a) is a constant depending only on s,n,a. 

This estimate combined with the assumptions p2.2p and /?£ > put into the equation (|12.4p 
gives that solves in the viscosity sense, and hence in view of the linearity also in the sense of 
distributions: 

< -(-A)'^/\^(x) < max{L(V', •)} + C'(s,"-,o-)||V'||l°° in (12.12) 

B3 

Hence are uniformly bounded in H"^ . Furthermore, as in the proof of [431 Proposition 2.9] 
(or alternatively using the Green's Function of (— A)*^/^ for ^3, see [M])) the boundedness of 
_(_A)'^/2^;£ also gives a uniform estimate for u^. Hence we can extract two subsequences 
u*" — > w, both locally uniformly and in H'^ weakly for some w. 

What remains to show is to confirm that w = v. We will show that w necessarily solves (|12.4p . 
and hence be uniqueness, w = v. We first note that since /?£ > 0, and hence by stability of 
viscosity solutions, also w, is a subsolution of F{w, x) > 0. Furthermore, /3 is chosen specifically 
so that F{i{j,x) < /3e(0) and hence is a supersolution of (|12.5p . Thus by comparison with , 

< ijj and hence w < ip as well. Therefore, w is a subsolution of (|12.4p . Now to justify the 
supersolution property suppose that w — (j) has a strict global max at xq G -B3. Then there exist 
Xe £ S3 with Xe xq and — (f) has a local max at Xg. If it happens that w{xq) = ^(xq), 
then (|12.4p is satisfied, and so we assume that w{xq) < iP{xq). In this case for e small enough, 
v^{xe) < V'Ca^e), and hence for e small enough, /3£(V'(xe) — u^(xe)) = 0. Thus F((/), Xg) < by the 
supersolution property of v^. Passing to the limit using the continuity of F and the smoothness 
of (j) gives F{(P,xq) < 0. Hence w solves (|12.4p . This concludes the proof of the lemma. □ 

To conclude this section, we give the brief proof of Proposition 112.11 



Proof of Proposition \12.1\ The inequality (jl2.12p in the sense of viscosity solutions and in sense 
of distributions is stable with regards to respectively local uniform limits and weak limits. 
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Therefore by Lemma Il2.5t v is in H"^ and hence for a.e.x, 

/ I .n+a dy < OO, 

JR" \y\ 

and V satisfies the inequality almost everywhere, and this gives the proposition. 



□ 
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